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Abstract 

Two-dimensional materials were first isolated no longer than ten years ago, and a comprehensive understanding of their properties 
under non-planar shapes is still being developed. Strictly speaking, the theoretical study of the properties of graphene and other 
two-dimensional materials is the most complete for planar structures and for structures with small deformations from planarity. The 
opposite limit of large deformations is yet to be studied comprehensively but that limit is extremely relevant because it determines 
material properties near the point of failure. We are exploring uses for discrete differential geometry within the context of graphene 
and other two-dimensional materials, and these concepts appear promising in linking materials properties to shape regardless of 
how large a given material deformation is. A brief account of additional contributions arising from our group to two-dimensional 
materials that include graphene, stanene and phosphorene is provided towards the end of this manuscript. 
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Geometry is a pillar of Science, and many physical theo¬ 
ries are decidedly geometric m This manuscript provides an 
overview of recent developments towards linking the properties 
of two-dimensional materials to a given arbitrary shape, where 
shape is understood as the local two-dimensional geometry of 
atom-thin materials that are embedded on a three-dimensional 
space. The unifying point of the description concerns the in¬ 
troduction of a discrete geometry to deal with two-dimensional 
materials while fully preserving their atomistic information. 

Thus, we showcase here a set of geometrical principles that 
apply to nets, where a net is a discrete surface or a mesh. We 
identify two-dimensional materials with meshes, and apply pre¬ 
cepts from a branch of Mathematics O that deals with discrete 
surfaces. We have presented a number of results in linking this 
geometry to materials properties already EH Eli 13. This 
subject has contributions from other teams as well ISlIglfTOl. 

I acknowledge Drs. Pacheco SanJuan, Wang, Harriss, Rivero, 
Vanevic, and Terrones’ contributions to this ongoing work. I 
also thank students Sloan, Horvath, Utt, Pour-Imani, Mehboudi, 
and Klee for their contributions at different stages. I am grateful 
to many colleagues for their observations and encouragement; 
most particularly to Maria A. H. Vozmediano. 

Graphene and other 2-D materials provide a stage to fur¬ 
ther our understanding of Physics. Perhaps the most natural 
connection to be studied concerns the creation of gauge fields 
on effective Dirac particles in 2-rl dimensions as the geome¬ 
try evolves from a reference, planar shape EiEiiia, to be 
addressed next. 

The starting point for us was the analysis of strain created 
by a scanning tunneling microscope (STM) on graphene m. 
There is a theory laid out on a structural continuum 1221 [TTI 
Ha [mill El [H [13 that correlates structural deformations to 
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mechanically-induced gauges on Dirac fermions in 2-rl dimen¬ 
sions. These effective Dirac fermions arise from a first nearest 
neighbor tight-binding description of 7r-electrons on graphene 
at low energies. Changes in distances arising from a structural 
deformation are estimated from a continuum model of the dis¬ 
tortion, and these changes in distances alter the magnitude of 
the tight-binding hopping terms locally. 

The formulation is inherently semi-classical, in the sense 
that the underlying dynamics is that of pseudospins (which strictly 
speaking are only valid on the ideal non-deformed crystalline 
structure) and the gauge fields produced by mechanical defor¬ 
mations induce local modifications to the 7r-electron pseudospin 
Hamiltonian. We estimated gauge fields employing that formal¬ 
ism (m[121 [131 but this question quickly came up: 

1. An STM can tell individual atoms. Can one rewrite the 
theory expressed on a continuum structure to reflect such 
atomistic nature? What do we learn when the theory is 
laid out this way that is different from the continuum for¬ 
malism? 

This paper contains three sections that are somehow inde¬ 
pendent: (1) Its main thrust is the description of the coupling 
of finite displacements to a semiclassical pseudospin dynam¬ 
ics of TT-electrons on graphene in which we attempt to provide 
answers to the two questions above (pages 2-5). (2) We then 
provide a description of a discrete geometry that applies to arbi¬ 
trary two-dimensional materials (pages 5-7). (3) The document 
ends by briefly mentioning other developments in graphene and 
other materials in which we have been involved (pages 8-10). 

1. A lattice gauge field theory for Dirac fermions in graphene 

The interplay among the electronic and mechanical prop¬ 
erties of graphene membranes remains under experimental and 
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theoretical investigation (TTJ [1210 [HI [20l [21 EH 1221 , and an 
insightful picture of the effects of deformations employs gauge 
fields that influence the dynamics of charge carriers 

EBEa. 

The formulation is inherently semi-classical and takes pseu¬ 
dospin hamiltonians as the main object, which strictly speaking 
are only valid on the ideal non-deformed crystalline structure, 
with gauge fields arising from slow-varying mechanical defor¬ 
mations providing local modifications to the said Hamiltonian. 

But graphene can sustain elastic deformations as large as 
20% m and using this picture, the resulting pseudo-magnetic 
fields are much larger than those magnetic fields available in 
state-of-the-art experimental facilities. The presence of a pseudo- 
magnetic field is observed via broad Landau levels (LLs) in 
strained graphene nanobubbles on a metal substrate (25\ . In 
addition to the pseudo-magnetic vector potential A^, strain also 
induces a scalar deformation potential Eg EB |26l EH that af¬ 
fects the electron dynamics in non-trivial ways. Part of our 
motivation was to reconcile the experimental results that can 
be obtained when the lattice is largely deformed, with a the¬ 
ory that by construction applies to small deformations. What 
we accomplished is a close view at the inner workings of this 
theory that has led to unique insights, and a quantitative under¬ 
standing of “slowly varying deformations” within the context of 
this theory. Our formulation brings to the spotlight some of the 
inherent assumptions on the prevailing theoretical framework. 

The underlying assumptions of the theory expressed on a 
structural continuum are expressed in the following sentence: 
“If a mechanical strain varies smoothly on the scale of inter¬ 
atomic distances, it does not break sublattice symmetry but rather 
deforms the Brillouin zone in such a way that the Dirac cones 
located in graphene at points K and K' are shifted in opposite 
directions m.” 

Previous statement tells us that -provided strain preserves 
sublattice symmetry- one can understand the effects of mechan¬ 
ical strain on the electronic structure in terms of a semiclassi- 
cal approach, in which mechanical strain induces the spatially- 
varying gauge fields Bs{y) = VxA^(r) and ^’^(r) into a spatially- 
varying pseudospin Hamiltonian Elps{(\, r), where is the 

low-energy expansion of the Hamiltonian in reciprocal space in 
the absence of strain. This semiclassical approximation is jus¬ 
tified when the strain extends over many unit cells and it pre¬ 
serves sublattice symmetry |[T2l[T3l[22l, and many of the equa¬ 
tions on this Section will help us keep track of said sublattice 
symmetry. 

Evidently, it is possible to determine the electronic proper¬ 
ties directly from a tight-binding Hamiltonian El in real space, 
without resorting to the semiclassical approximation and with¬ 
out imposing a sublattice symmetry a priori. That is, while the 
semiclassical Elpsi^i, r) is defined in reciprocal space (thus as¬ 
suming some reasonable preservation of crystalline order), the 
tight-binding Hamiltonian 7T in real space is more general and 
can be used for membranes with arbitrary spatial distribution 
and magnitude of the strain. 

In the previous formulation of the theory both and Es are 
expressed in terms of a continuous displacement field u(v,y) 
obtained within first-order continuum elasticity (CE) ETJ [121 



Figure 1: (a) Definitions of geometrical parameters in a unit cell, (b) Sublattice 
symmetry relates to how pairs of nearest-neighbor vectors (either in thick, or 
dashed lines) are modified due to strain. These vectors change by Ary and 
At' upon strain (7 = 1,2). Relative displacements of neighboring atoms lead 
to modified lattice vectors; the choice of renormalized lattice vectors will be 
unique only to the extent to which sublattice symmetry is preserved: At'. ^ At^. 


[T3l|22l- It is not possible to assess sublattice symmetry on a 
continuum media, and therefore proper phase conjugation of 
pseudospin Hamiltonians becomes an implicit assumption of 
that theory. 

But the only way to know whether the strain preserves sub¬ 
lattice symmetry |[T2l implies analyzing relative atomic displace¬ 
ments in arbitrary structural distortions that could be captured 
directly from experiment ll28ll , or from molecular dynamics sim¬ 
ulations. 

Let us start by considering the unit cell before (Eig. ^sl)) 
and after arbitrary strain has been applied (Eig.[2b)). The lat¬ 
tice vectors and the vectors joining atoms are given by (Eig.[2a)): 

ai = (1/2, V3/2)ao, 82 = (-1/2, V3/2)ao, (1) 

Vs’ 

( 2 ) 

before the deformation takes place. When a deformation is ap¬ 
plied (Eig.[2b)) the two off-diagonal terms making up the pseu¬ 
dospin tight-binding Hamiltonian will be: 

3 

j=i 


J=i 

where t is the hopping term, St is its change upon strain to be 
explicitly defined later on, and k is the crystal momentum. 

Each local pseudospin Hamiltonian will only have physical 
meaning when it is properly conjugated, which implies sublat¬ 
tice symmetry holds. This happens at unit cells where: 

At;. - At7 0=1,2). (3) 
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One also notes that on this first nearest neighbor picture the 
diagonal terms are always real even when At'. 9 ^ Atj (more on 
this later). 

Condition ^ can be re-expressed in terms of changes of 
angles Aaj or lengths ALj for pairs of nearest-neighbor vectors 
Tj and t'j [j = 1 is shown in thick solid and 7 = 2 in thin dashed 
lines in Fig.[ 2 b)]: 

{Tj + Atj) • (Tj + At') = \tj + AtjWtj + At)| cos(AQfy), (4) 

sgn(AQfy) = sgn ([(t; + Atj) X (tj + At')] • k ), (5) 

where k is a unit vector along the z-axis, sgn is the sign function 
(sgn(a) = - 1-1 if (2 > 0 and sgn(a) = -1 if a < 0 ), and: 

ALj = \rj + ATyl - \Tj + At)|. (6) 

Previous expressions indicate that the sublattice symmetry 
Ifni does not hold a priori. In the continuum approach, both 
ATj and At', are captured at the same point in space using an 
identical value of the deformation field u(v,y) and hence the 
structural aspect just uncovered is hidden. Forcing this symme¬ 
try to hold in the lattice depicted at Figure lb amounts to im¬ 
posing an artificial mechanical constraint 1291 , and we re-derive 
the theory without using that continuum deformation field. The 
details of the discrete model we developed follow. 

In the absence of mechanical strain, the reciprocal lattice 
vectors bi and b 2 are related to the lattice vectors by |[30l : 

(7) 


where ^ and B = with the 

\<221 < 222 / \b 2 l ^ 22 / 

choice we made for ai and a 2 : 


1 \ 2 ;: 
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bi = 1, — —, and b 2 = -1, — — 


1 \ 2;7r 
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( 8 ) 


The ^-points on the first Brillouin zone are given by (c.f., 
Fig.|3a)): 


K, = K2 


^ and K 3 


bj -i- 2 b 2 


(9) 


and: 

K 4 = -Ki, K 5 = -K 2 , and K 6 = -K 3 . (10) 

The relative positions between atoms change when strain is 
applied: Tj Tj At^ (j = 1,2,3), and -Tj -Tj - At'. 
( 7 = 1 , 2 ). 

How do reciprocal lattice vectors change under mechanical 
load (to first order)? Taking Eqn.j^as the starting point, Aaj and 
ALj must all be close to zero. In that case we set At'. ^ At j 
for j=1,2 and continue our program. We then define: 

Aai = Ati - At 3 , and Aa 2 = At 2 - AT 3 , ( 11 ) 


or in terms of (two-dimensional) components: 


AJl = 


(Atu - At^x 

\ATiy - At^j 


At2x - AT3;,\ 

AT2y - ATsy) 


( 12 ) 


(a) No strain: (b) Under mechanical 

strain: 



Figure 2: First Brillouin zone (a) before and (b) after mechanical strain is ap¬ 
plied. The reciprocal lattice vectors are shown, as well as the changes of the 
high-symmetry points at the corners of the Brillouin zone. Note that indepen¬ 
dent K points (K and K') move in the opposite directions. The dashed hexagon 
in (b) represents the boundary of the first Brillouin zone in the absence of strain. 


The matrix Jl changes to A2R, and we must modify B 

so that Eqn. 0 still holds under mechanical load. To first order 
in displacements becomes: 

+ ^ (13) 


By comparing Eqns. (7) and ( [T^ , the reciprocal lattice vectors 
in Fig-I^b) must then be renormalized by: 

AS = . (14) 


This additional term is evident when working directly on the 
atomic lattice but it was missed in Ref. lED. Let us now calcu¬ 
late shifts of the ^-points due to strain. For example, K 2 (= K 
in Fig.j^a)) is shifted by: 


AK = AK 2 = 


3«o 
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Ati;,: + At 


2 x 


■2At 


3jc 


2X9 


V3 


and using Eqn. (10) one gets AK' - -AK 2 , so that the K (K 2 ) 
and K' (-K 2 ) points shift in opposite directions, as expected 

caEii. 


1.1. Gauge fields induced by strain 


Equation ^ tells whether mechanical strain varies smoothly 
over interatomic distances ina. This observation provides the 
rationale for expressing the gauge fields without ever leaving 
the atomic lattice: When At'. ^ At, at each unit cell a mechan- 

j 

ical distortion can be considered “long-range,” and the first- 
order theory is valid. Local gauge fields can be computed as 
low energy approximations to the following 2 x 2 pseudospin 
Hamiltonian: 



(15) 


withg = + dO)e'(’’^+^’’A(K»+AK„+q)^ ^ ^ 

Keeping exponents to first order we have: 


(Tj -f ATj) • (K„ -f AK„ + q) - Tj • K„ -h Tj • AK„ -h At^ • -f Tj • q. 

The exponent is next expressed to first-order on At: 

^i(Tj-Kn +Tj-AKn +ATj-Kn +Tj -Q) ^ 

je'T-K,, T^. . q + e'T-K,, [ I + + Atj • K„)] . (16) 
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Carrying out explicit calculations one sees that: 

3 

y] + iijj ■ AK„ + At,- • K„)] =0. (17) 

For example, at ^ = K 2 we have: 


1 + 


Ai7l{^Tu + At2x + AT3;,) 
9ao 


(1+^3 -^ 3 ), 


with phasors adding up to zero. Similar phasor cancelations 
occur at every other ^-point. 

The term linear on AK„ on Eqn.[T^cancels out the fictitious 
^-point dependent gauge fields proposed in Ref. 1^ . which 
originated from the term linear on Atj on this same equation. 
Equation ( p3] ) takes the following form to first order at K 2 in 
the low-energy regime: 


"Hps = 


+ 


( 0 


( 18 ) 


with the first term on the right-hand side reducing to the stan¬ 
dard pseudospin Hamiltonian in the absence of strain. The 
change of the hopping parameter t is related to the variation 
of length 01221: 


6tj = — -Tj ■ At,-. (19) 

% 

This way Eqn. ( p^ becomes: 

"Kps = hvpfr • q + (), ( 20 ) 

with /* = ^[2 t3 • At 3 - Ti • Ati - T 2 • At 2 + V3/(T2 • At 2 - 

Ti • Ati)], and hvp = 2^^. The parameter /i can be expressed 
in terms of a vector potential: As f\ = -hvp^. This way: 


1/7100 r 2T3-AT3-Ti-ATi-T2-AT2 

7:al *- V3 

-/(T2 • At2 - Ti • Ati)], 


( 21 ) 


The deformation potential can be written in terms of the 
average (Edef) and the difference (Emass) between Eg, a and Es^b 
(Eqns. ( [22| ) and ^23\ ) at any given unit cell: 

Edef — 2 ^Es,a + Es^b\ and Efnass ~ \{Es,a-Es,b\ (24) 

Both quantities are of the order of tens of me Vs. It is worth 
noting that many other teams do not include the deformation 
potential in their models, even though it can lead to significant 
changes in the electronic spectrum ElEI. 

We also note that, while the determination of required 
the preservation of sublattice symmetry, the diagonal terms Es^a 
and Es^b are not constrained by this requirement because the 
diagonal entries in Eqns. 22 and 23 remain real regardless of 
the actual magnitudes of At j and of At'.. 

Equations (21-23) are discrete gauge fields', this is, they take 
a single value at any given unit cell. These Equations thus 
provide an original, discrete, viewpoint -in which one retains 
atomic positions- to the issue of gauge fields in two-dimensional 
materials Eiiiiia. 

The mass term, Eqn. ( |24l i leads to a Zeeman term that arises 
as a second-order difference relation among potential energies 
for an atom on the A-sublattice at the ^-point (c.f. Fig. [DEI 

El [33 El: 

^ ((<5rf - drf>) - {6tf - dtf) (25) 

+ (drf-drf)-(<-<)). 

Here, /ig is the Bohr magneton (^ 5.8 x 10“^ eV/Tesla), 6tJ^ 
is the standard change in hopping upon strain at unit cell n = 
1,2,3 iniiniiii, and ^2.5. The pseudomagnetic field 

ffle 

Bs changes sign at the B-sublattice and/or at the K' point. Eg 
is the average deformation potential at a given unit cell (H aris¬ 
ing from the rearrangement of the electron cloud upon strain 
I 22 I . Equation ( [25] ) provides a “microscopic” vehicle to obtain 
the local magnitude of the pseudo-magnetic field at any given 
unit cell directly. We refer readers to Publications EHma for 
explicit calculations of gauge fields and electronic spectra in 
graphene membranes with specific shapes. 


and the diagonal entriesin Eqn. are deformation po¬ 
tentials that arise on a neutral system as it is deformed. They 
indicate the change of the local electronic density as the system 
is distorted: 


o.3^y 1 


|Tj- - AtjI - ap/ V3 
ao/V3 


( 22 ) 


and 


Es,b - - 


03eV 1 y I'^'i - - «o/ V3 

0.12 3 E ao/V 3 


(23) 


The deformation potential as expressed above has been taken 
to linear order on the average bond increase following explicit 
results from ab-initio calculations EH. 


7.2. Relation to the continuum formalism 


The continuum limit is achieved when ^ 0 (for / = 

1,2, 3). We have then (Cauchy-Bom rule): TyArj ) 

where Uij are the entries of the strain tensor, and Eqn. ( j^ be¬ 
comes: 


^100 

2^/3na{ 


■ ( ^xx ^yy 2iUxy), 


(26) 


-J’ 


as expected lfT^[T3]| . 

Edef is an average over changes of distances, and hence re¬ 
flects the basic form found in terms of the deformation tensor 
El Edef ^ Uxx + Uyy, while E^ass takes its continuum form in 
Eqn. ( |25] ) |[3^ : reference lE) contains further details. 

We have thus addressed the first question raised in Page 1: 
We have rewritten the original theory taking into consideration 
the deformation at individual unit cells. 
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Figure 3: (a) The finite-difference curl leading to the pseudomagnetic field Bg 
[Eq. is obtained from hoppings among an atom on the A-sublattice and 
three neighboring atoms on B-sublattices. 


One of the exciting points of graphene is that it furnishes a 
field theory for Dirac electrons in 2+1 dimensions. The theory 
we work with involves fermions having a (pseudo-)spin arising 
from the 7r-electrons on two inequivalent sites on a honeycomb 
lattice. Unlike an intrinsic spin, a unit cell can be thought as a 
plaquette that has a finite spatial extent. Each of these plaque- 
ttes can be assigned two integer indexes (/, j) that hence furnish 
a discrete lattice. On a crystal with N atoms, there are N/2 such 
unit cells that hence define a discrete space. As we work on 
directly on this discrete space, we say that we realize a discrete 
field theory. The next point to cover is the discrete geometry in 
which this lattice gauge theory takes place. 


2. The discrete geometry of two-dimensional materials 


Similar to the statements made in previous section, graphene’s 
geometry is commonly studied in terms of the continuous dis¬ 
placement field u = as well. Specifically, on thin- 

plate continuum elasticity the strain tensor is = (daUp + 
dpUa + daUydpUy + daZdpz)l2, with z an out-of-plane elongation 

There, differen¬ 
tial geometry and mechanics couple as: 


gap - ^ap + 2Uap, 


Kp - n 




(27) 


where is a tangent vector field, dap is the reference 

(fiat) metric and h = is the local normal |[T3l|23l|26l[32l 

[33l . However, peculiarities of how graphene ripples 1381 [^l40l 
smu, slides and adheres (391 na may be beyond first-order 
continuum elasticity. 

We have investigated alternative geometrical frameworks to 
deal with discrete atomistic surfaces. This is an important en¬ 
deavor because geometry is behind the spin diffusion in rippled 
graphene ||44l [45l, behind the chemical properties of confor¬ 
mal (non-planar) two-dimensional crystals 0 , and may even 
herald the strain engineering of two-dimensional crystals with 
atomistic defects, an area completely unexplored so far. 
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Figure 4: Schematics of the parameters employed to determine the local dis¬ 
crete geometry for graphene. Ap is the small hexagon colored in marine blue 
representing the area of a (Voronoi) unit cell, and the angles 0/ add up to 2n on 
a fiat local surface. 


Such discrete geometry exists EEllTl. There, the Wigner- 
Seitz/Voronoi unit cells that span a locally-evolving area Ap are 
the underlying discrete geometrical objects, and the atomistic 
information is always preserved. The discrete formalism for ge¬ 
ometry rests on interatomic distances without a mediating con¬ 
tinuum, just as the theory for the electronic response of Dirac 
fermions did in previous section. In what follows, we review 
the tools for geometrical analysis and study the local geometry 
of rippled graphene 1461 . 

The four invariants that determine a local shape arise from 
the metric (g) and curvature (k) tensors as follows E|6l: 

Tr(g), Det(g), H = Ti:(k)/2Tr(g), K = Det(^:)/2Det(g), (28) 

where Tr (Det) stands for the trace (determinant), H is the mean 
curvature and K is the Gaussian curvature, respectively. 

The discrete metric is defined from the local lattice vectors 
a„ EEl gap = ^a ' ^p^ and the discrete Gauss curvature (K^) 
originates from the angle defect l5l [6l [9l l47l l48l : 

6 

KD = (27T-J]0dlAp. (29) 

i=l 

Here 6i (i = 1, ...,6) are angles between vertices and Ap will 
be defined below; see Fig. The Voronoi tessellation gener¬ 
alizes the Wigner-Seitz unit cell to conformal two-dimensional 
geometries ElISl. 

The discrete mean curvature Hj) measures relative orienta¬ 
tions of edges and normal vectors along a closed path: 

6 

Hd = Y^ e,' X (v,;,+i - Vi-ij) ■ h/4Ap. (30) 

i=l 

Here, V/ is the position of atom i on sublattice A, and C/ = V/ - 
is the edge between points p and i (note that ai( 2 ) = ei( 2 )). Vij+i 
is the normal to edges e, and e^+i (i is a cyclic index), and n = 
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Figure 5: (Color online) (a) Creation of ripples by cutting a square with side 
L = 0.21 iJ.m at 1 Kelvin: The membrane trades a planar configuration for a 
rippled one. (b) Geometrical invariants within the dashed square shown in (a). 

is the area-weighted normal with Af = |e/ x e/+i |/2 121. 

2ji=l 

and Ap = ^ Yj%i 

The discrete metric and curvatures furnish geometry con¬ 
sistent with a crystalline structure and they lead to the faithful 
characterization of graphene’s morphology beyond the effective- 
continuum paradigm, Eqn. ( [27] ). This is advantageous when the 
atomic conformation is known from molecular dynamics (e.g, 
ll46l ) or experiment (e.g., 1^ ) because: (i) fitting of the atomic 
lattice to an effective continuum is not needed any more, (ii) 
the Chemistry of conformal graphene can be addressed from 
the discrete geometry El and, since atoms are always avail¬ 
able, (iii) the discrete theory brings new insights and under¬ 
standing into the physical theory (e.g., non-preservation of sub¬ 
lattice symmetry, the form of gauge fields O, the creation of 
mass from strain EiEa). We emphasize that the discrete ge¬ 
ometry is accurate regardless of elastic regime, hence it can be 
used to verify whether the conditions for continuum elasticity 
hold in the problem at hand. 

2. 1 . The geometry of rippled graphene 

The importance of a sound geometrical framework is moti¬ 
vated by rippled graphene. We contrast ripples created by ther¬ 
mal fluctuations ll46ll with those created at low temperature due 
to edges. These two mechanisms lead to different types of ge¬ 
ometries (hence different magnitudes of strain-derived gauges). 
In a system with periodic boundary conditions, thermal fluctu¬ 
ations create significant changes in interatomic distances (i.e., 
in the metric) ll^ and -as the boundaries are fixed- such in¬ 
creases on interatomic distances produce out-of-plane deforma¬ 
tions (i.e., rippling). 

Now consider a square graphene sample with about three 
million atoms, in which strain was relieved at the low tem¬ 
perature of 1 Kelvin. The resulting membrane is shown in 
Fig. 5(a), where colors indicate varying heights across the sam¬ 
ple (H. Ripples in Ref. ||46]| originate from increases in the 
metric. On the other hand, the white margin in between the 
“rippled” (curved) sample and the (yellow) exterior frame high¬ 
lights an apparent contraction of our finite sample when seen 
from above. 


The details of this geometry are shown in Fig. 5(b): Det(g) 
and Tr(g) are unity almost everywhere (yet there are significant 
random fluctuations driving the scales). The metric tells us that 
the membrane does not contract, and its area thus remains al¬ 
most unchanged. We show in Fig. 5(b) the discrete curvatures, 
highlighting cusps by ovals, valleys by squares, and ridges by 
triangles. Cusps and valleys have the largest Gaussian curvature 
Kb (deep red), while ridges have the smallest one (deep blue). 
As expected, the mean curvature Hj) takes its largest (smallest) 
value at valleys (cusps) and alternates sign around ridges. The 
curvature -without metric increases- explains the white mar¬ 
gins on Fig. 5(a). 

The discrete geometry reflects the mechanism leading to 
ripple formation, thus highlighting the virtue of a geometry that 
originates from atoms. An accurate determination of Hj) is im¬ 
portant since Hu leads to spin diffusion in rippled graphene 
1441 l45t . Though much has been said about ripples, no geo¬ 
metrical study with the accuracy provided in Ref. Elia exists. 

The starting point in the continuum theory is a fiat metric 
Sap. There, a non-zero curvature directly leads to increases in 
interatomic distances Eqn. and a non-zero height is di¬ 
rectly identified with a non-zero strain-derived gauge. A ques¬ 
tion then arises whether the sample under study actually obeys 
Eqn. ^21) . The situation shown in Fig. 5 is a counterexample 
to the geometry inferred from Eq. •Hz), because the metric is 
almost constant, even though the height profile z is clearly non¬ 
fiat. Gauge fields for similar samples were reported in Ref. IH . 
Fig. 5 represents the accurate geometrical characterization of 
rippled graphene down to the atomic level. We studied the ge¬ 
ometry of graphene under load, providing lattice gauge fields 
and electronic spectra, in Ref. l|6|. 

We next address another aspect of conformal two-dimensional 
materials: The potential increase in chemical reactivity under 
non-planar, conformal shapes. 

2.2. Chemical measures and geometry 

We first provide two known chemical measures for carbon- 
based materials: 

1. Mean Bond Fength (MBF): Aromaticity is not a directly 
measurable property and hence it cannot be defined un¬ 
ambiguously. Yet, the structural representation accom¬ 
modating the maximum number of Clar sextets best rep¬ 
resents chemical and physical properties, and Clar sextet 
migration increases chemical reactivity. How aromatic is 
conformal graphene l[46l ? Can rippling be explained in 
terms of the creation of “aromatic domains”? The mean 
bond length (MBF) is defined as follows l[5l l49l[50ll : 

MBL = dec - 7 X! (31) 

where acc,i are bond lengths on a closed loop El. Ac¬ 
cording to Ref. EqI, MBL is a reliable tool for analysis 
of large aromatic systems. 
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(a) MBL (Ang) (b) Tr(g) 



1 1.394 
1.392 
1.390 
1.388 
1.386 



1 0.968 
0.966 
0.964 
0.962 
0.960 
0.958 


(c) Det(g) 


1 0.934 
0.930 
0.926 

1 0.922 
0.918 




Figure 6: Chemical measures and the local geometry of a rippled two-dimensional crystal. While no direct correlation can be drawn among MBL in (a) and metric 
measures in (b-c), the piramidalization apyr in (d) is directly proportional to the mean curvature H in (e). Actual atomic bonds are seen in black O. 


2. The angle Og-n between the bonds and the normal vector 
n at atom p has a single value Ocm under a spherical ge¬ 
ometry ED and Ocrn = 7112 on di flat surface. 6cm can be 
generalized for arbitrary geometries as an average: 

1 ^ 

0crn = TYj^{i)crn, (32) 

i=\ 

with 6{i)o-T: the angle among a bond vector and the local 
normal. Equation ( [3^ takes its usual form for fullerenes, 
where 0{i)cm = Og-n for all bonds ED- The pyramidaliza- 
tion angle apyr was deflned by Haddon as follows ED: 

CKpyr = Kn “ :^/2. (33) 

The degree of sensitivity of MBL with respect to fluctua¬ 
tions on interatomic distances makes a direct correlation diffi¬ 
cult Il46l [6ll4l. In Figures 5(a-c) we contrast MBL with Tr{g) 
and det{g) (we will not display the determinant of the metric 
tensor in further figures). Details of the creation of the rippled 
structure can be found in prior work E [H |6l . Recalling the 
notion that pristine graphene has equal bond lengths, Fig. 6(a) 
indicates that atomistic fluctuations will have a bearing on the 
aromatic behavior of rippled samples; this concept has not been 
discussed before, nor its ramifications. 

On the other hand, there exists a remarkably simple, one-to- 
one correlation among the pyramidalization angle and the mean 
curvature (sign included) for all the systems studied, as already 
evident from Fig. 6(d-e): 

apyr (in rads) ^ 1 x Hd (in A“^). (34) 


Equation ( [34| ) is an interesting result because it relates a com¬ 
monly used angular measure for orbital hybridization and chem¬ 
ical reactivity with the mean curvature. This result was hidden 
in plain sight; this shows once again how the discrete geometry 
makes plenty of sense. 

apyr is a signed quantity, as follows: Direct inspection of 
Eqn. ( [^ indicates that apyr will be positive for a bulge, and 
negative for a sag. Similarly, the mean curvature Ho -Eqn. ([30|)- 
is a vector quantity projected onto the local normal; the relative 
orientation of the normal (facing “up” or “down”) confers Hu 
with a sign as well. (Geometrically speaking, one sees that ra¬ 
dius of curvature changes sign for a bulg or a sag, so Hu must 
be signed.) But the correspondence goes beyond the sign. The 
cross products on H -Eqn. confers an additional sinu¬ 

soidal function, which approximates as the angle rather well up 
to 20 degrees (0.35 rad), within the range of all pyramidaliza¬ 
tion angles we saw. The correlation given by Eqn. is re¬ 
markable as it informs our intuition concerning hybridization, 
thus making the mean curvature a direct tool for analysis of 
hybridization and chemical reactivity for two-dimensional sys¬ 
tems with 5- and p electrons. 

Additional results in Ref. El include the discussion of a 
geometry of systems with topological defects such as fullerenes 
Eai EH mi ESI Ei HD m, Schwarzites 1591 . ionic crystals 
|60||6T1, and other hexagonal systems with atomistic defects. 

This concludes the main discussion on this contribution. As 
advances in two-dimensional materials are occurring at great 
speed, we end this work briefly highlighting some contributions 
to such endeavor arising from our group. 
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3. Brief account of additional contributions on graphene 

and on other two-dimensional materials 

The tools employed in previous sections are a combination 
of electronic structure within a tight-binding approach, molec¬ 
ular dynamics, and a discrete geometry. 

But the group has been invested on models of electron trans¬ 
port that capture the electronic structure of graphene and metal 
leads with a Green’s function approach that is coupled to an 
ab-initio electron Hamiltonian (^[631 [64l; these models pro¬ 
vide noise features that reproduce experimental features ll65l 
not seen in more basic theoretical models 16611. 

Another result from our group that is becoming relevant 
and is worth mentioning concerns the experimental observa¬ 
tion on an STM of a bulk material through graphene and its 
theoretical confirmation 1^; the result is interesting because 
it can be used, for instance, to study the surface of black phos¬ 
phorus through graphene or through hexagonal boron nitride 

monolayers. (Ml 

We briefiy discuss in what follows our contributions to other 
two-dimensional materials. The techniques employed in these 
studies combine ab-initio methods, tight-binding models, and 
the discrete geometry. 

3.1. Stanene 

Proceeding by direct analogy to silicene and germanene 16^ . 
known studies of the electronic properties of two-dimensional 
tin (7011211 [m have been performed under the implicit assump¬ 
tion that the HB phase is not viable. Contrary to this assump¬ 
tion, we determined using ab-initio methods that the HB two- 
dimensional structures of heavy column-IV elements tin and 
lead are stable and lower in energy than their LB counterparts 
(c.f., Fig.|^), thus representing the true optimal structures of 
these two-dimensional systems (731. The HB structure is a 
hexagonal close-packed bilayer (c.f.. Fig. [7^). 

Haldane’s honeycomb model has been studied in closed ge¬ 
ometries ca and one of the many candidates for its practical 
realization is LB tin (stanene). Unfortunately, a fullerene-like 
Sn6o is not stable (Fig. [7]:) so tin and lead are no-go elements 
for topological fullerenes ca 

We also determined ca that the optimal phase of two- 
dimensional fiuorinated stanene is not analogous to tetrahedrally- 
coordinated graphane l75l as it was postulated in Refs. (TTHtTI 
either. There is no indication for tetrahedral coordination of 
tin atoms in bulk fiuorinated tin (76l and tetrahedral coordi¬ 
nation (TT] |72l does not yield the most stable two-dimensional 
fiuorinated tin either. 

Indeed, as seen in Fig. [^, the phase space for decorated 
two-dimensional tin is larger than originally anticipated: The 
graphane-like phase (7TJ[72]| realizes the metastable minima la¬ 
beled 6 that turns into phase 4 upon in-plane compression. In 
the optimal structure, 7, four-fold coordinated Sn atoms form a 
sequence of parallel zig-zag one-dimensional chains with two 
fluorine atoms mediating interactions among neighboring Sn 
chains. The structure is realized on a triangular lattice with 
oq = 5.230 A [Fig.[^]. The Wigner-Seitz unit cell is within the 



Figure 7: (a) The high-buckled (HB) phase is more stable than the low-buckled 
(LB) phase, making freestanding stanene ^2^2) metastable, (b) The actually 
stable, high-buckled phase, is a HCP bilayer with trivial electronic properties (it 
is a metal), (c) As seen on this one-minute-long structural optimization, stanene 
does not realize topological fullerenes l74l either. 

dotted area in Fig. [^, where the symmetry axes are shown as 
well. 

The first Brillouin zone in Fig. [^a) shows a top view of 
the conduction band and the high-symmetry points in momen¬ 
tum space. As seen in Fig. [^b), the arrangement of parallel 
ID Sn wires gives rise to an electronic structure with only two 
anisotropic Dirac cones on the First Brillouin zone located away 
from the K-points at positions Vi and V 2 = +0.85Ki, respec¬ 
tively. From now on we identify the v-axis with the line join¬ 
ing tin atoms across fiuorine bridges. The Fermi velocity is 
close in magnitude to that of graphene and it is anisotropic: 
vpy = 5.4 X 10^ m/s [Fig. 9(c)], and = 2.1 x 10^ m/s 
[Fig. [^d)] and a 2A = 0.02 eV gap opens due to SOI, five 
times larger than the intrinsic gap due to SOI in graphene da. 
Phase 6 transitions from a topological insulator to a trivial in¬ 
sulator ca, but the electronic structure of the optimal phase 
remains robust under larger isotropic strain. 

The electronic dispersion in Fig. |^b-d) can be understood 
in terms of a 2 x 2 7r-electron tight-binding Hamiltonian (77]| 
in which an effective coupling f is set among the tin atoms 
originally linked by fiuorine bridges [thin bonds on Fig.j^b)], 
and t is the coupling among actual Sn-Sn atoms [thick bonds on 
Fig. [^b)]. Using interatomic distances among Sn atoms from 
Table I we obtain the blue dashed lines in Fig. |^c,d) with t = 
0.8 eV and t' = —t which reproduce first-principles results. 

To account for SOI, we realize an oblate low-energy Dirac- 
Hamiltonian at the vicinity of the Vi ,2 points. The numerical 
results on Fig.j^e) are consistent with a coupling r^o-^Sx (73]| . 
Indeed, eigenvectors of TzCTxSx project spins onto the -v, -fx, 
-fx, -X axis parallel to the Sn-F bonds, inverting signs at each 
valley and lacking sublattice polarization, consistently with ab- 
initio data [Fig. 9(e)]. Thus, the low-energy dynamics is given 
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Figure 8: (a) Phases of two-dimensional fluorinated tin; structures shown to 
the right, (b) Symmetries of the most stable structure (7), depicting triangular 
(dashed) and Wigner-Seitz (within dotted perimeter) unit cells, the two symme¬ 
try axes, and the two Sn sublattices A and B. Structural stability is demonstrated 
by (c) phonon dispersion curves and (d) the structural stabilization of a finite- 
size sample. 



by: 

H = -ih^^VpxTzO-A + VFyO-ydy)'^ + 'V\AT:,crxSx)^. 

An unprecedented specific coupling of momentum -including 
direction- with spin oriented along x and valley degrees of free¬ 
dom is thus realized by the second term in previous equation. 
The valley degree of freedom can be addressed by a bias along 
the Vi - V 2 axis that breaks inversion symmetry. Similarly, 
a magnetic field along the x axis will break time-reversal sym¬ 
metry, locking the valley and crystal momentum direction at the 
Vi, V 2 points. The dynamics invites the use of two-dimensional 
fluorinated tin for valleytronic applications. 

The structural stability of HB tin and HB lead clearly have 
fundamental consequences for the practical realization of substrate- 
free non-trivial topological phases based from these elements. 

3.2. Phosphorene 

Studies of planar phosphorene with defects have begun to 
appear |[78l [791 and single-digit-percent strain on planar black 
phosphorene induces a ten-fold change on its electronic gap 
|[80l . Curvature can induce strain ED, and we indicated how 
to induce a positive Gaussian curvature on black phosphorene, 
how to characterize such geometry, and how this shape infiu- 
ences the electronic gap. Additionally, we demonstrate that this 
reduction in the fundamental gap applies to other phosphorene 
allotropes. 

Finite-size phosphorene cones -with a positive Gaussian 
curvature- were built by removing the angular segments seen 
in Fig. 10, joining atoms along the red lines, and a subsequent 


Figure 9: (a) Conduction band on the first Brillouin zone, highlighting high- 
symmetry points and locations of valleys Vi and V 2 away from the K- and 
K’-points, (b) The two valleys on the Brillouin zone arise from the two-fold 
symmetry of the atomic structure, (c-d) Band structures along high-symmetry 
lines, including a two-band tight-binding fit. (e) Spin texture resolved over 
valley (t), energy, and sublattice (cr) degrees of freedom. (The spin projection 
onto the z-axis is of the order of 1% at most.) 

atomistic relaxation with molecular dynamics at the ab-initio 
level. (An area segment must be added instead to build struc¬ 
tures with negative Gaussian curvature -e.g. ; Elia.) The an¬ 
gular sections removed subtend a 46° angle for black phospho¬ 
rene. To get a suitable joining line for black phosphorene, the 
starting point is a planar structure with a dislocation at an an¬ 
gle of 26.8° that does not create localized electronic states lIT^ 
but confers additional structural rigidity. The edges were passi¬ 
vated with Hydrogen atoms, and the conical structure on Fig. 10 
contains about five hundred atoms. The electronic gap for the 
planar fiake takes a constant value of 1.1 eV, as highlighted by 
the yellow color in the last column of Fig. 10. The magnitude 
of the gap is about twice its nominal value (0.4 eV in standard 
density-functional theory liSTlfO ) due to finite-size effects. 

The black phosphorene cone subtends a solid angle of 0.83x27r 
radians, thus inducing a positive Gaussian curvature, and it shows 
compressive strain at the apex as indicated by the white tones 
in their metric around this point. Black phosphorene has a large 
structural rigidity due to its ridged structure and as a result the 
metric and the curvatures in conical structures lack a perfect ra¬ 
dial symmetry. The discrete metric and Hj) tell us the locations 
of the dislocation line in the planar structure. 
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Figure 10: The discrete geometry and the semiconducting gap of black phos- 
phorene cones. The cones are built by removing the indicated segments on an 
initially planar structure. A remarkable ~20 percent reduction of the semicon¬ 
ducting gap occurs due to the strain induced by the conical shape. 

The last column in Fig. 10 depicts the semiconducting gap 
at each atomic position, and the darkest color indicates a 20% 
reduction of the gap with respect to its value in the planar struc¬ 
ture, due to the curvature-induced structural compression (SD 
discussed in previous paragraph. Thus, topological defects can 
help in tuning the local gap of phosphorene Q. 

4. Conclusion 

In conclusion, we have presented here a number of con¬ 
tributions to two-dimensional materials in which their discrete 
geometry plays a preponderant role. 
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